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New  branching  rules  and  new  methods  to  control  the  size  of  LP's  resulted  in  a 
world's  record  for  the  solution  of  large  traveling  salesman  problems. 


15  NUMBER  OP  PAGES 


16  PRICE  CODE 


17.  SECURITY  CLASSIFICATION  18  SECURITY  CLASSIFICATION  19  SECURITY  CLASSIFICATION  20  LIMITATION  OF 
OF  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UNCLASSIFIED  SAR{SAME  AS  REPORT) 


98  3  4  02  2 


StJ^'cta'd  'o-— 

•  t5^  *v,  ..  .  -7 

;  -fm  ■ ' :  2 


7540'0' -280*5500 


Final  Report  on  AFOSR  Grant  ^ F496A0~92~J-0053 


Investigations  in  Linear  and  Integer  Prograrnining 
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Linear-Programming  Tools  in  Integer  Programming: 
The  Traveling  Salesman  Problem 


A  new  "world  reconl  "  w,(s  ,rt  ui  ih«*  Uavdirig  salfMuan  prfibietn  I  lus  work, 
joint  with  D.  Apples'atc.  \  ('hvalal  and  W .  (’cMik.  w'as  designated  h_v  Dis 
cover  Magazine  as  on**  <d’  it>  tup  tesi  >cience  stories  of  the  year,  i  he  AFOSR 
support  of  work  on  lineai  •i>rograinining  tools  for  integer  programming  pro¬ 
vided  partial  support  for  this  ISP  project. 

The  problem  we  solveii  had  dOfS  nodes.  I'he  total  computation  time  wa.s 
equivalent  to  appfo.vimalely  one  and  one-half  years  of  SPARC  2  seconds.  The 
previous  “record”  wa.s  2d92  nodes,  a  problem  that  we  can  now  solve  in  under 
two  hours  on  a  SPARC’  2  (With  no  branching!}.  In  addition  to  solving  the 
3038  problem,  vve  solved  15  additional  previously  unsolved  real-world  TSP 
instances  (mostly  VLSI  drilling  problems). 

VVe  are  now  in  the  process  of  preparing  to  write  up  the  results  of  the 
above  work,  a  project  that  we  hope  will  be  completed  within  the  next  six 
months.  From  the  point  of  view  of  integer  programming,  the  impact  of  this 
work  should  go  well  beyond  'he  LSP.  Several  new  separation  techniques  were 
developed  along  with  new  met  hods  for  controlling  the  sizes  of  the  the  linear 
programs  that  result  in  am.  t>ranch-and-cut  approach.  Other  consequences 
were  improved  general  pm  i-us.>  !. ranching  rules  and  improvements  in  the  in¬ 
teraction  of  the  LP-soh  ei  a  u  h  t  lie  branching  code  (not  to  mention  an  almost 


two-orders-of-magni  tilde 
in  this  class). 


!,jm>veinent  in  the  LP  solutions  times  for  problems 
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Basis  Recovery 

A  joint  paper  wilii  Matt  Sail/tiutn.  "Ut-i  uvt-ri!ii'  an  i>pli!nai  I. P  b.iM-  froii!  an 
interior  point  solutiosi.  "  a.i-,  .  .Uui  subunlltni  for  puhh<  atjoii  In 

paper  we  consnier  a  pinn.il  u!  tiie  ba»i.^  rti-ovtry  prcjiflcin  for 

point  Uiethods;  (lir  en  an  Mnf  nn.i!  j>!  nna!  MtlntiOn  x*  fur  Mime  hin'af  jirograiii. 
construct  an  o[.>tiinai  fni'in  /f*  I  tu'  essentials  of  that  proi'etiure  aia^, 

•  (,’onstruct  some  '  nutre-t  i<  .isiv  stable”  starting  b;iais  ti  tliat  iin  ludes  its 
many  "positise”’  e'riu.-iits  liom  s'  a,s  possible. 

•  After  the  consti  tterjun  .4  tlie  initial  B.  there  will,  in  general,  remain 
some  “positive"  (ijinponents  of  s'  that  are  outside  of  B.  Such  variables 
are  designated  as  .>upf  rba^u:.  I'lie  secotid  part  of  the  procedure  is  to 
modify  the  simplex  method  to  deal  with  these  superbasic  tariables. 

Results  from  this  approach  v%-ere  far  superior  to  those  for  any  previous 
basis- recovery  procetiuriv  In  tests  on  the  .N'etlih  problem  s<.*t.  this  method 
was  able  to  recover  an  optimal  Irasis  in  5'a,  of  the  interior-point  .solution  time 
for  most  problems. 

Since  the  completiim  ot  this  work  there  have  been  further  significant  de¬ 
velopments.  In  particular,  a  procedure  proposed  by  Megiddo  has  been  inves¬ 
tigated  and  shown  to  i>e  very  promising  for  problems  that  are  highly  primal 
degenerate.  This  proceilun'  ha.s  a  time-cornple.xity  that  is  linear  in  the  total 
number  of  rows  and  column'  m  the  given  linear  program. 


Basis  Reduction  for  Integer  Programming 

The  object  of  this  resea r-  h  wa.'  to  study  a  method  that  is  particularly  well 
suited  for  mixed  integer  luogramming  problems  with  a  relatively  small  num¬ 
ber  of  general  (non  0/1  i  mtet/cr  variables.  The  drawback  of  this  method  has 
been  that  the  computatiuu.d  cifort  grows  quickly  with  the  number  of  integer 
variables  in  the  problem  i.cm^  solved.  Thus,  with  current  technology,  the 
number  of  integer  varial>I«  s  is  essentially  restricted  to  100  or  fewer.  The  aim 
of  our  research  has  been  to  extend  t  he  applicability  of  this  method  by  study¬ 
ing  how  it  can  be  applied  in  parallel.  In  the  process,  we  have  not  simply  been 
aiming  to  implement  a  known  algorithm  in  a  new  computing  environment. 


but  to  investigate  the  cfiaracler  of  tfu?  ineliuui  when  (ojHuirent  step',  ate 
executed. 

The  method  m  questiun  wa.s  developed  m  ih'JO  by  Cmik.  Kuthrrfurd. 
Scarf  and  Shallcross,  io  dale,  we  have  been  able  to  devekjp  a  »et|uent!al 
implementation  of  their  a!y,urithm  that  is  vapabie  uf  seh. ing  ^ojne  iiuali  but 
very  difficult  mixed  integ.-r  psogratiis  that  ario-  in  tfie  design  uf  fiber  uptn 
networks.  These  probiems  .tre  nit i actable  wuli  standard  firanch  and  Bound 
approaches.  Our  sequentia!  (  vale  has  b«'<ni  tested  un  tlie  Intel  Istjij  fi>  j»etrube. 
and  the  preliminary  steps  needed  to  ubtain  a  parallel  implementation  are 
nearly  complete. 

The  algorithm  that  we  have  impleiiientesi  procwds  by  creating  a  large 
search  tree,  and  using  a  very  complicated  routine  to  determine  which  branch 
to  follow.  Our  focus  thus  far  has  been  on  parallelizing  the  branching  step  of 
the  overall  process,  and  it  is  our  belief  tlial  this  will  drastically  reduce  the 
time  required  for  solving  a  problem.  Otice  this  step  has  been  completed,  we 
plan  to  investigate  possible  parallelization  of  the  complicated  routine  that 
controls  the  branching  sequence.  This  routine  is  invoked  within  every  node 
of  the  large  search  tree  that  is  generated  by  the  algorithm,  and  we  hope  that 
by  finding  a  viable  parallelization  of  this  routine,  the  running  time  of  our 
code  will  be  reduced  enough  that  wc  will  be  able  to  solve  problems  with  a 
larger  number  of  integ^^r  variables. 

Polyhedral  Approaches 

This  work  focuses  on  the  study  of  the  underlying  polyhedral  structure  and  the 
development  of  a  branch-and-cut  IP  solver  for  a  class  of  structured  0/1  integer 
programs  arising  from  a  truck  dispatching  scheduling  problem.  The  problem 
is  characterized  by  a  group  of  set  partitioning  constraints  and  a  group  of 
knapsack  equality  constraints  of  a  specific  form.  Families  of  facets  for  the 
polytopes  associated  with  iiuiividual  knapsack  constraints  are  identified,  and 
in  some  cases,  a  complete  ( haracterization  of  a  polytope  is  obtained.  In 
addition,  a  notion  of  ■■conflict  graph"  is  introduced  and  utilized  to  obtain  an 
approximating  node-packing  poiytope  for  the  convex  hull  of  all  0/1  solutions. 
We  are  in  the  process  of  developing  the  branch-a’^d-cut  solver  which  generates 
cuts  based  on  both  the  knapsack  constraints  and  the  approximating  node¬ 
packing  polytope,  and  incorporates  these  cuts  into  a  tree-search  algorithm 
that  uses  problem  reformulation  and  LP-based  heuristics  at  each  node  in 
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the  search  tree  to  assist  in  tlie  soiuluiii  tnocess.  Nunierical  exi>erinients  will 
then  be  perfornit^ci  on  large-scale  real  instances  su|j|)lied  by  an  Oi!  (unipanv 
in  the  United  States.  1  he  sjjitinia!  schedules  obtained  torrespoinl  to  cost 
savings  for  the  troiupativ  and  gre.tter  job  satisfaction  for  drivi..,  dste  to  more 
balanced  work  schedules  and  income  distribution,  it  is  noteworthy  that  this  is 
apparently  the  first  lime  that  branch  am!  cut  has  bwn  applied  to  an  ec|ualily 
constrained  problem  in  which  the  entries  in  the  constraint  matrix  and  right 
hand  side  are  not  purely  0  1 
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